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Abstract 

We formulate a covariant transport approach for high energy nucleus-nucleus 
collisions where the real part of the nucleon selfenergies is htted to nuclear matter 
properties which are evaluated on the basis of a NJL-type Lagrangian for the quark 
degrees of freedom. The parameters of the quark-model Lagrangian are hxed by the 
Gell-Mann, Oakes and Renner relation, the pion-nucleon S-term, the nucleon energy 
as well as the nuclear binding energy at saturation density po- We hnd the resulting 
scalar and vector selfenergies for nucleons to be well in line with either Dirac- 
Brueckner computations for p < 2po or those from the phenomenological optical 
potential when accounting for a swelling of the nucleon at hnite nuclear matter 
density. The meson-baryon interaction density is modelled to describe a decrease of 
the meson mass with baryon density. The imaginary part of the hadron selfenergies 
is determined by a string fragmentation model which accounts for the in-medium 
mass of hadrons in line with the ’chiral’ dynamics employed. The applicability of 
the transport approach is demonstrated in comparison with experimental data from 
SIS to SPS energies. The enhancement of the K“’“/7r“’“ ratio in A -|- A collisions 
compared to p -|- A reactions at AGS energies is reproduced within the ’chiral’ 
dynamics. Furthermore, detailed predictions for the stopping in Pb -|- Pb collisions 
at 153 GeV/A are presented. 


1 Introduction 


The study of hot and dense nuclear matter by means of relativistic nucleus- 
nucleus collisions is the major aim of high energy heavy-ion physics. However, 
any conclusions about the nuclear properties at high temperature or baryon 
densities must rely on the comparison of experimental data with theoretical 
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approaches based on nonequilibrium kinetic theory. Among these, the covari¬ 
ant RBUU approach [1-9], the QMD [10] or RQMD model [11] have been 
successfully used in the past. As a genuine feature of transport theories there 
are two essential ingredients: i.e. the baryon (and meson) scalar and vector 
selfenergies - which are neglected in a couple of approaches - as well as in¬ 
medium elastic and inelastic cross sections for all hadrons involved. Whereas 
in the low-energy regime these ’transport coefficients’ can be calculated in the 
Dirac-Brueckner approach starting from the bare nucleon-nucleon interaction 
[12,13], this is no longer possible at high baryon density (pB A 2-3po) and high 
temperature, since the number of independent hadronic degrees of freedom in¬ 
creases drastically and the nuclear system is expected to enter a phase where 
chiral symmetry is restored [14-17]. Such a phase transition is dynamically 
due to a change of the nonperturbative QCD vacuum at high temperature 
or baryon density and the chiral invariance of the interaction between quarks 
and gluons in the QCD Lagrangian. As a consequence the hadron selfenergies 
in the nuclear medium should change substantially especially close to the chi¬ 
ral phase transition and any transport theoretical study should include the 
generic properties of QCD that so far are known from nonperturbative com¬ 
putations on the lattice [18-21]. However, such nonperturbative calculations 
will not be possible for high baryon densities within the next years and we 
have to rely on suitable effective Lagrangians that lead to the same physical 
condensates and thermodynamic behaviour as the original QCD problem. 

In this paper we aim at formulating a ’chiral’ transport theory for the hadronic 
degrees of freedom, which in covariant notation formally can be written as a 
coupled set of transport equations for the phase-space distributions /h(a;,p) of 
hadron h [1-3,7,8], i.e. 


{(n, - n.ay-," - a;} Mx.p) 

= f (72(73(74 . . . [Gr^Gr]i2^34..Ar(n + n 2 — Hs — n 4 ...) 

h2h3h4...'^ 

X {fhsix,P3)fw{x,P4)fhix,p)fh^ix,p2) 

- fh{x,p)fhAx,P2)fhs{x,P3)fhAx,P4)} ■ ■ ■ • (1) 

In eq. (1) U^{x,p) and U0{x,p) denote the real part of the scalar and vector 
hadron selfenergies, respectively, while [G'“’“G'] 12 ^ 34 ...(^r(n + 02 — n 3 — 04 . . .) 
is the ’transition rate’ for the process 1 + 2 —)• 3 + 4 + ... which is taken 
to be on-shell in the semiclassical limit adopted ^ . The hadron quasi-particle 


^ The index F at the ^-function indicates that off-shell transitions of width F should 
also be allowed. In the actual transport simulation, however, we use the on-shell limit 


F = 0. 
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properties in (1) are defined via the mass-shell constraint [7], 




( 2 ) 


with effective masses and momenta given by 


Af(((a;,p) = Afh + Ul{x,p) 


Ili^{x,p)=p'^ -Ul^{x,p) , 

(3) 

while the phase-space factors 


fh{x,p) = 1 ± fh{x,p) 

(4) 


are responsible for fermion Pauli-blocking or Bose enhancement, respectively, 
depending on the type of hadron in the hnal/initial channel. The dots in 
eq. (1) stand for further contributions to the collision term with more than 
two hadrons in the hnal/initial channels. The transport approach (1) is fully 
specihed by U^{x^p) and U{^{x,p) (p = 0,1,2,3), which determine the mean- 
held propagation of the hadrons, and by the transition rates G^GS‘^{. . .) in the 
collision term, that describe the scattering and hadron production/absorption 
rates. 

The scalar and vector mean helds and f//) are conventionally determined in 
the mean-held limit from an effective hadronic Lagrangian density £h which is 
the sum of the Lagrangian density for the free helds and some interaction 
density i.e. 

rH = E£f. + £“ ■ (5) 

h 

The actual form of however, is only known for more simple cases at 

low baryon density pB and its general form at high pB and for large relative 
momenta between the interacting hadrons - which is probed in nucleus-nucleus 
collisions up to 200 GeV/A - is essentially undetermined. This opens up a large 
parameter space for coupling constants phh', from factors at the vertices as well 
as respective powers in the hadron helds, which might lead to various density 
isomers in the nuclear equation of state or mesonic condensates, respectively. 

In order to reduce this large parameter space and to incorporate aspects of 
chiral symmetry we here adopt the strategy to specify for baryons via 
an effective Lagrangian for the underlying quark degrees of freedom £q for 
nuclear-matter phase-space conhgurations (see below). In hxing scalar and 
vector couplings as well as vertex cutoffs for the ’quarks’ and by comparing 
the energy density for nuclear-matter conhgurations from £q with that of £h 
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on the hadronic side, we can ht hadronic couplings and vertices in even for 
high baryon densities and thus determine and in the transport equation 
(1) in a less arbitrary way. 


The ’hard’ hadronic processes, on the other hand, which govern the r.h.s. of 
eq. (1), are modeled by the LUND string-fragmentation [22] which is known 
to describe inelastic hadronic reactions in a wide energy regime at zero baryon 
density. The medium modihcations due to the hadron selfenergies, however, 
require to introduce some conserving approximations in line with the density 
dependent hadron masses. With the specihcations of U^{x^p) and U{^{x,p) and 
the inelastic collision rates G^GS{...) the transport approach (1), which will be 
denoted by HSD ^ , is fully dehned and can be confronted with experiment. 


Our work thus is organized as follows: In Section 2 we will evaluate the nucleon 
selfenergies basis of a NJL-type Lagrangian Cq and model the 

meson-baryon interaction. In this respect we hrst hx the free parameters of 
our quark model interaction by the Gell-Mann, Oakes and Renner relation 
[23], the pion-nucleon S-term and the free ’nucleon’ mass. We then extend 
our study to the computation of the energy density of nuclear matter conhg- 
urations, discuss the necessary modihcation of the nucleon formfactor in the 
medium and present results for the nuclear equation of state at low and high 
baryon density. Scalar and vector selfenergies for the nucleons are obtained by 
htting the coupling parameters in the covariant approach of Weber et ah [7] 
to the density dependence of the effective mass and the energy per nucleon 
from the NJL calculations. We compare our results to the low density, low 
energy Dirac phenomenology and discuss the extrapolations to high baryon 
density and large relative momenta. Furthermore, the meson-baryon interac¬ 
tions in are described along the line of Kaplan and Nelson [24] essentially 
employing density-dependent meson masses. In Section 3 we specify the modi¬ 
hcations of the familiar LUND string-fragmentation model [22], which models 
the imaginary part of the hadron selfenergies or collision rates, to include (as 
a hrst step) the modihcation of the hadron masses at hnite baryon density. In 
Section 4 we apply our transport approach to nucleus-nucleus collisions from 
SIS to SPS energies and test its applicability in comparison with experimental 
data. As a hrst test for the partial restoration of chiral symmetry in heavy- 
ion collisions we compute the ratio for systems at AGS energies and 

compare to the available data. Section 5, hnally, is devoted to a summary and 
discussion of open problems. 


^ Hadron-String-Dynamics 
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2 Hadron selfenergies 


In this section we specify the evaluation of the mean helds and that 
enter the l.h.s. of the transport equation (!) for the mean-held propagation. 
In order to reduce the parameter space and to obtain extrapolations for 
and at high baryon density, we use an effective Lagrangian density Cq for 
quarks that is compatible with the approximate chiral invariance of the QCD 
Lagrangian. We thus hrst hx the effective Lagrangian Cq for the quark degrees 
of freedom on the mean-held (one loop) level for low energy QCD problems, 
where the gluon helds A“(a;) are supposed to be integrated out. The effective 
interaction determined in this way should not be used in further perturbation 
theory (e.g. for scattering or transition rates) since it is assumed to be the 
result of an inhnite resummation of interaction diagrams. In a second step we 
then extract the real part of nucleon selfenergies from quark conhgurations 
that describe nuclear matter at hnite density. 


2.1 The effective quark Lagrangian 


The underlying idea of an effective 4-point interaction for quarks has already 
been discussed e.g. by Vogl and Weise in ref. [25]. Since the fundamental 
currents in QCD are color currents, i.e. J“ = an elementary color 

current interaction with a universal coupling Gc is expected to be dominant. 
An effective Lagrangian for 8{xi — a; 2 )-like quark interactions thus reads 

- niQffq -GlY^ , ( 6 ) 

a = l 


where C(a = 1,...,8) are the 5'17(3)coior matrices with = Sab/^^rho a 

diagonal mass matrix in havor space, i.e. rho = diag{m^, m[|, m°) and tpq = 
(u,d,s) is the quark spinor in case of 5'17(3)fiavor- The color-current interac¬ 
tion is invariant under chiral transformations or 5'17(3)fiavor rotations. The 
Lagrangian (6), however, in its present form is not yet well suited for the for¬ 
mulation of quark dynamics on the mean-held level because antisymmetriza- 
tion generates a further mixing of color, havor and Dirac indices. It is thus 
more convenient to introduce a Fierz transformation, i.e. to antisymmetrize 
the 4-point interaction to proceed with further computations on the Hartree 
level. The Fierz transform then generates color singlet as well as color octet 
terms, i.e. [25,26] 


Cq{x) = tljq{iffdf, - mo)V’q + 


8 = 0 



+ 


V’qihsyyV’q 
, 2 ! 
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2 


-^v I] I 

-G'c I] (^q7A*^“V’q)" - ‘^Gl ^ ^ I [^qY^“V’q] + ^^q^S Y^“ V’q^ ^ 

a=l a=li=0 y \ / \ / j 

+ ^^C I (^^q7A*Y^“V’q^ + (^^q7A*75Y^“V’q^ | , (7) 

where Gg = 2Gy = |Gq. In (7) the matrices A*(z = 1,..,8) stand for the 
*hG(3)fiavor degrees of freedom with fr(AjAj) = 26ij while A° is given by 
A° = with I 3 denoting the 3x3 unitary matrix in flavor space The 

Lagrangian (7) in its color-singlet version has been the starting point for RPA- 
type calculations for the bosonic excitations of the nonperturbative QCD vac¬ 
uum, i.e. the mesonic degrees of freedom [25-27]. Similar Lagrangian densities 
have also been exploited by a variety of authors [28-36] following an early 
suggestion by Nambu and Jona-Lasinio (NJL) [37]. 

In our present study we will discard the mesonic (RPA-type) sector and con¬ 
centrate on the determination of a static effective quark-quark interaction by 
nucleon properties as well as nuclear matter related quantities. Similar con¬ 
cepts have been proposed by Guichon [38] and Saito and Thomas [39] based 
on bag-model wavefunctions. Here we start with a slightly different concept 
by determining the quark wavefunctions for the nucleon from the experimen¬ 
tal data for the proton electromagnetic formfactor. In this way we intend to 
circumvent the problem of absolute conhnement which cannot be dealt with 
properly using only a color neutral mean-held approach of the NJL-type. Since 
we are only interested in energy densities for given quark conhgurations, the 
resulting Lagrangian should not be used for dynamical studies such as the 
RPA response (mesonic sector). Furthermore, it is not expected that the re¬ 
spective soliton solution of (7) for a nucleon presents a dynamically stable 
object of a shape consistent with the experimental proton formfactor. 


2.2 The isospin symmetric nucleon system 


In this Subsection we concentrate on vacuum as well as nucleon properties, 
where the nucleons are assumed to be represented by 3 additional valence 
quarks with a hxed phase-space distribution on top of the (truncated) Dirac 
sea with a formfactor in line with the experimental data. The singlet terms of 

^ The most general four-point interaction compatible with QCD symmetries starts 
from combinations of all possible vector and axial currents. Therefore, in general, 
there is no strict relationship between Gs and Gy. 
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the Lagrangian (7) in the mean-held limit - performing the sum over the havor 
matrix elements - then leads to the following Lagrangian for = (u,d,s), 


- ^k) V’fc 

A;=u,d,s 

f - 2 - 2'] 

+^|(V’fcV’fc) + (V’fcWsV’fc) I 

+ (^fc757^V’fc) }}, (8) 

where the couplings and Gy are now considered as free parameters. For the 
systems of positive parity, which are of interest in our present work, also the 
pseudoscalar and pseudovector terms vanish in the Hartree limit such that we 
are left with the scalar and vector term, only. This is quite similar to the cr-io 
model [40,41] in the nuclear physics context. The hamiltonian density then is 
given by 


T~('{x) = {V’fc (-17*^* + "1° - G's(V’fcV’fc)) V’fc 

A;=u,d,s 

H ^ 5 (9) 

which leads to the gap equations for the effective masses m/;, i.e. 

rrik = ml-Gl {tpk'tpk)- ( 10 ) 

Since the problem (9) decouples in the flavor degrees of freedom we will con¬ 
sider in the following only u-quarks assuming = m[| and neglect a possible 
strangeness content of the nucleon furtheron. 

For the nonperturbative vacuum we then end up with the gap equation in 
phase space for the effective quark mass of u or d quarks: 

= m° + GIjG^ I d^p 0(As- | P |) = my, (If) 

(27r) J /p2 _|_ ^2 


where we have introduced a spatial cutoff parameter As to regularize the 
divergent integral over the Dirac sea. Alternatively, one might also introduce 
covariant cutoff schemes as in [25,34], but for reasons to be discussed below 
in context of eq. (13) we prefer to use the scheme (11), since we are basically 
interested in quark conhgurations with a well dehned rest frame. In eq. (11) 
the factor g = 6 arizes from the trace over color and spin in eq. (10). The 
gap equation (11) then leads to a constituent quark mass in the 

nonperturbative vacuum. 
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The coupling constant Gs together with the cutoff parameter As now can 
be determined via the Gell-Mann, Oakes and Renner relation [23] assuming 
(uu) = (dd), 

K fl = -« + (uu) , (12) 


where = 93.3 MeV is the pion decay constant, the physical pion mass 
and (uu) the scalar condensate (for u or d quarks in the vakuum). Choosing 
= 7 MeV as an average value of the light quark mass the quark condensate 
then amounts to (uu)^/^ sa —230 MeV; a value which is achieved by choosing 
a cutoff As ~ 0.59 GeV and Gs ~ 4.95 GeV“^ in (11). 

In the presense of additional localized light valence quarks on top of the Dirac 
sea the gap equation (11) modihes locally to 


mu(r) = m° - Gg 
.2 9 


3 u^u(r) 


d^p 


+Gi 


(27r)3 J ^p2 _|_ y72^(r)2 

mu(r) 


/u(r,p) 


(2 


TT 


d'^p 


Jp^ + mu(r)2 


0(As- I P I), 


(13) 


where /u(r, p) denotes the phase-space distribution of a single u-quark which 
has to determined in a model dependent way. 


In a fully dynamical theory on the mean-held level /u(r, p) should result from 
the solution of the Dirac equation 

- G'sPs(r) + 7°pv(r)}V’fc(r) = 7°^fcV’fc(r) (14) 


with 


ps(r) = (V’fc(r)V’fc(r)) , 

dv(r) = (V’I(r)V’fc(r)) , (15) 

and subsequent Wigner-transformation of V’l-('”~s/2)V’fc('” + s/2)- However, 
since we do not aim at a dynamical theory for the nucleon - due to the lack 
of conhnement in Cf^{x) (8) - and we are only interested in the total energy of 
well dehned quark conhgurations, we hx /u(r, p) (from outside) by the experi¬ 
mental electromagnetic formfactor of the proton which is well represented by a 
dipole approximation up to momentum transfers 25 GeVyG [42]. This 

implies that the quark charge distribution (of a proton) is of the exponential 
form [43] 

(V’q(r)V’q(r)) sa iVo exp{- \ r | /bo) = Pq(r), (16) 
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where r is given in fm, bo = 0.25 fm and A^o = provides normalization 

to 1. Considering now a nucleon state averaged over spin and isospin, i.e. a 
mixture of proton, neutron and A's of average mass Mn ~ 1.085 GeV, we 
obtain for the u-quark density 

3 

du(r) ^ -No exp{- \ r | /bo), (17) 

where the factor 3/2 reflects the average u-quark content of the states con¬ 
sidered. In the local density approximation the phase-space distribution for 
u-quarks then is given by 

/u(r,p) = 0(pF(r)- I p I) (18) 

with the local Fermi momentum 

pp(r) = (Q/g 7r2)^/Vu(r)^/h (19) 


This approximation has been quite successfully applied in the nuclear physics 
context [1,2] and also been adopted in [44,45] for quark oriented models. It is 
a legitimate approximation for the quark phase-space distribution as long as 
one is interested in expectation values like the total energy, only. 


Inserting /u(r,p) (18) with (17) and (19) in the gap equation (13) we can 
compute the effective quark mass mu(r) for the ’nucleon’ described above. The 
resulting coordinate-space dependence of mu(r) (full line) for the ’nucleon’ is 
shown in Fig. 1 together with the u-quark density (u^(r)u(r)) = Pu{r) (dashed 
line). In the interior of the ’nucleon’ the effective quark mass drops to about 
= 7 MeV and thus the quark scalar selfenergy f/g to zero. 


Whereas the scalar sector now is hxed by the gap equation (13) for arbitrary 
quark phase-space distributions /u(r, p) - that are at rest within the frame of 
reference considered here - the local vector quark interaction is modihed in 
order to allow for an explicite momentum dependence. We note that nonlocal 
generalizations of the NJL Lagrangian have been suggested by Bowler and 
Birse [46]. We adopt a similar concept and assume that the vector interaction 
in (8) is mediated by massive color neutral (vector) gluons which implies to 
modify the couplings 


Gy 


Gy 


Ay + 


( 20 ) 


where Ay sa 1.2 — 1.5 GeV is a vector cutoff and q denotes the momentum 
transfer in the quark-quark interaction. This strategy is similar to that used 
in effective meson-exchange interactions for hadron-hadron scattering [47]. 
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Fig. 1. Effective quark mass m = m-a(r) (full line), quark density p-air) (dashed line) 
and scalar condensate — < qq (dotted line) as a function of the radial distance 
r from the center of the ’nucleon’. 

The energy density r™(r) in phase-space representation thus reads (including 
a factor of 2 from the summation over u and d quarks) 

r°°(r )=25 j /^(r, p) 

-2^ J + "^u(r)2 0(As- I p I) 

+2 |tGi^s(.-)^ + j‘ep.<ep2k(v,v,) ^ p^)2 /»(--.P2)} 

-^vac, (21) 

where the vacuum contribution 

= -2 h/ 0 (As- I p I) + (22) 

has been subtracted. In eq. (21) ps = {^-a — fn^)IG\ is the scalar quark density 
and the vector quark density is pv = ( 2 ^^ / /u ('”5 p)- The total energy {H) 

of a quark conhguration described by /u(r, p) then is obtained by integrating 
/ d^r r°°(r). 



10 







The average nucleon energy to be fixed in our case corresponds to 1.085 GeV, 
which is the average of the nucleon and the A mass. Since in eq. (21) for 
r°°(r) all quantities are determined except the quark vector coupling Gy and 
cutoff Ay, the vector coupling (for fixed Ay sa 1.5 GeV) is well determined by 
the total energy of the quark configuration. Our fit provides Gy = 4.2 GeV“^ 
using /u(r, p) = 0(pF(r)— | p |) with PF(r) from (19). The pion-nucleon 
S-term, defined by the following matrix element with the nucleon state, 

SvrN = + ?7^d) (N|uu + dd|N) , (23) 

within the parameters stated above leads to S^n ~ 47 MeV, which is well in 
line with the value extracted from pion-nucleon s-wave scattering of 45 ± 7 
MeV from [48]. This will be of significant importance for the scalar nucleon 
selfenergy later on. 

We stress again that the coupling parameters Gs, Gy and cutoffs As, Ay only 
apply for the semiclassical static quark configurations discussed so far and 
should not be considered as appropriate for a fully dynamical theory on the 
basis of the Lagrangian (8). In fact, the vector coupling Gy is larger than in 
refs. [25-27] where the mesonic sector has been explored. As a consequence 
we can only attempt to describe ’nucleons’ at finite baryon density and have 
to discard mesonic degrees of freedom. 


2.3 Symmetric nuclear matter 


In order to evaluate the energy density for symmetric nuclear matter configu¬ 
rations we have to introduce in addition to /q = /u(r, p) a phase-space distri¬ 
bution for the nucleons or ’localized’ quark states /n. Denoting by (fn, Pn) the 
position and momentum of a nucleon, the corresponding quark phase-space 
distribution /q(r, p)rN,PN obtained from a translation of the center of /q by 
Fn and a proper Lorentz transformation by /3 n = in phase space, 

i.e. a contraction of /q by = ^1-/3^ in coordinate space and dilation 
in momentum space by 7 n, which keeps the individual phase space integral 
invariant. 

Before going over to the nuclear matter problem we first consider two-nucleon 
configurations for ’frozen’ nucleon quark distributions in the nucleon-nucleon 
center-of-mass system (c.m.s.). As an example the local quark phase-space dis¬ 
tribution /u(f, p) - as met in the overlap regime of two colliding quark states - 
is depicted in Fig. 2 as a function of and for py = 0. It’s macroscopic pa- 

This is denoted as the ’sudden’ approximation in the nuclear physics context. 
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Fig. 2. Characteristic quark phase-space distribution in the overlap regime for two 
colliding nucleons for = 0. The sphere with radius As characterizes the Dirac sea 
contribution at rest. 

rameters are given by the relative momenta Pi , P 2 of the ’quark wave functions’ 
with respect to the nuclear matter rest frame and the individual Fermi mo¬ 
menta ppi (r), PF 2 (r) that are determined by the individual densities at space 
position r by PFi(r) = {6 /as before. The Dirac sea contribution 
at rest is indicated by the sphere with radius As. 

Due to the 3-momentum cutoff As in the gap equation (13) our prescription 
is not Lorentz-invariant and all quantities computed depend on the reference 
frame. For nucleon-nucleon collisions the natural frame of reference is the 
c.m.s., i.e. Pi = P 2 , which should be small compared to the nucleon mass. 
For the following illustration we thus restrict to small relative momenta of the 
’nucleons’ P = 3(Pi + P 2 ) = 6 P 1 . 

Iteration of the gap equation (13) for the two nucleon system then yields 
the effective mass mu(r; Pi, Pi, p^) as well as the scalar density ps(r)- The 
resulting quark vector (solid lines) and scalar densities (dashed lines) are dis- 
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played in Fig. 3 as a function of 2 ; = r for x = y = 0 for different distances R 
of the two nucleons (and a constant relative momentum P = 0.2 GeV/c). Due 



r [fm] 


Fig. 3. Spatial quark distribution 2/3 Pu{x = 0,y = 0,2; = r) (full lines) and scalar 
quark density (dashed lines) for two colliding nucleons with relative momentum 
P = 0.2 GeV/c for different relative distances R from 0.8 fm to 2.4 fm. 

to the gap equation (13) the scalar quark density drops substantially even for 
a moderate overlap of the nucleons [R sa f.Ofm), which reflects the ’interme¬ 
diate’ range (cr-held) attraction of the two nucleons, whereas the overlap of 
the vector densities becomes more substantial at short distance [R sa 0.5fm), 
which reflects the cu-held in terms of the conventional cr-io model [40]. 
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We note that due to the non-covariant cutoff As in the gap equation (13) 
the effective mass of a quark becomes momentum dependent for nucleon- 
nucleon conhgurations even in the c.m.s. This is more quantitatively shown 
in Fig. 4 where P) is displayed as a function oi Pi = P 2 = P and 



P 


Fig. 4. Effective quark mass m = P) as a function of P = Pi = P 2 and 

p^ = p^ = according to the gap equation (13). 

p^ = p^ = Pm to illustrate the smooth general dependence on density and 
relative momentum. It is clearly seen from Fig. 4 that the effective mass drops 
with density p and increases for hxed p with the relative momentum P^. = 2 P. 
Since the origin of this momentum dependence is not of dynamical nature - 
e.g. a hnite range of the scalar quark-quark interaction - one has to worry 
about its consequences for the nuclear-matter computations we aim at. For 
Pb = lOpo we get a nuclear Fermi momentum pp sa 0.57 GeV/c and thus 
for the relative momentum parameter P we have P = pp/3 ~ 0.19 GeV/c. 
A closer look at Fig. 4 then tells us that the effective quark mass is nearly 
independent on P for P < 0.2 GeV/c such that the momentum dependence 
of mu is rather insignihcant for our purposes. 

However, before evaluating the energy density for nuclear matter conhgura¬ 
tions we have to make sure that for vanishing nuclear density the energy of a 
nucleon moving with momentum pN = 3pu agrees with the dispersion relation 
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of a free nucleon, i.e. 


E{p^) = ^JpI + M2 , (24) 

where Mn is the nucleon mass in its rest frame. This is indeed the case as 
shown in Fig. 5, where the relation (24) (dashed line) is compared to the 
result from integrating r™(r) over r (solid line). The slight deviations from 



Pn [GeV/c] 

Fig. 5. Free nucleon dispersion relation (24) (dashed line) and / d?r r°°(r) (solid 
line) for a ’nucleon’ moving with momentum Pn. 

the exact result (24) provide a measure for the violation of covariance in the 
model adopted here which, however, are not serious for nucleon momenta 
Pn < 0.57 GeV/c. The result in Fig. 5 comes about because the total energy 
of the nucleon is dominated by the relative kinetic energy of the quarks. 

We now continue with the isospin symmetric nuclear matter problem where 
the nucleon phase-space distribution for hxed spin and isospin is given by 

/N(rN,PN) = 0 (pf- I Pn I) , (25) 

with the nucleon Fermi momentumpp = (O/dvr^pN)^^^, where pN is the nuclear 
matter density which will be discussed in units of po ~ 0.17 fm“^. 

Here a further problem is related with the change of the nucleon formfactor 
in the medium. As suggested e.g. by the interpretation of the FMC effect by 
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Close et al. [49] or arguments based on chiral symmetry by Brown [16] the 
nucleon might change its size in the nuclear medium such that the vector 
density of a quark is no longer given by (16). A fully dynamical model of the 
nucleon in the nuclear medium should give this modihcation of the formfactor 
in a selfconsistent manner. Since the Lagrangian (8) here is only considered 
to provide an effective quark-quark interaction for the energy density (21) we 
model such in-medium effects by modifying the width parameter bo in (16) as 

&o(/On) = 0.25fm ^1 + a —^ (26) 


with a parameter a to be determined by the nuclear matter saturation point 
(see below). 

In order to carry out computations for the nuclear matter problem we simulate 
the quark phase-space distribution /u(r,p) - which enters T™ in (21) - by 
characteristic samples k 

/u(r,p) = ^/u(r-CN,P-PN;^o(/?Jv)) , (27) 

j=i 


where /u(r, p, bo) denotes the semiclassical quark phase-space distribution for 
a ’nucleon’ of width 6 o(/On) (26). The nucleon positions are determined 
by Monte Carlo in a box of volume V = with —a/2 < < a/2. 

Only those samples are accepted for which the average distance to the next 
neighbour agrees within 3% with that for the respective inhnite nuclear matter 
value. The nucleon momenta then are selected by Monte Carlo with the 
constraint |p/^| < Pf{pn) and = 0- Additionally we rejects samples 

where the average kinetic energy 

Tn = 4 E ( ApTF+MS - Vn) (28) 


does not match with the nuclear matter value within 3%. The density pN in 
these simulations is given by pN = AjV (input) while A = 64 has been adopted 
throughout the calculations. In order to compile the dependence of the total 
energy on pN we have scaled the individual positions r/j with a ~ p^^^^ and 
the momenta p/j with pj/^. 

A snapshot of the quark density (for hxed z) for a chacteristic sample k at 
normal nuclear matter density po is shown in the upper part of Fig. 6; the 
resulting effective mass mu(r) according to the gap equation (13) - for the 
conhguration shown in the upper part of Fig. 6 - is displayed in its lower 
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part. Since at normal nuclear matter density the overlap of the nucleons is 
only moderate, the individual scalar ’quark bags’ can still approximately be 
separated in space for a given time. As an example for higher nucleon density 



X [fm] 



X [fm] 

Fig. 6. Snapshot of the spatial quark distribution (for hxed z) at normal nuclear 
matter density po (upper part) for a = 0.18 together with the resulting effective 
quark mass m = m-a(x,y) (lower part). 

we show a snapshot of the quark distribution at 4 X po for a = 0.18 in Fig. 7 
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(upper part) together with the corresponding quark mass m^{x,y,z =const) 
(lower part) from the gap equation (13). Since the overlap of the quark distri¬ 
butions now becomes substantial, the average quark mass drops to about 30 
MeV indicating partial chiral symmetry restoration. 



X [fm] 



X [fm] 


Fig. 7. Snapshot of the spatial quark distribution (for hxed z) at 4 X po (upper part) 
together with the resulting effective quark mass m = m-a(x,y) (lower part). 
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Now performing the integration of r°°(r) over coordinate space and averaging 
over characteristic samples k for nuclear matter conhgurations (as shown in 
Fig. 6) ® , dividing by the number of nucleons on the grid and subtracting the 
bare nucleon mass we can compute the energy per nucleon ~ 100), 


and thus establish a direct link between the energy density of quarks with the 
energy per nucleon of isospin symmetric nuclear matter at hnite density p/po- 
In (29) the energy density r™(r) is dehned by eq. (21) with /u(r,p) replaced 
by /u(r,p) from (27). 

The energy per nucleon (29) (for a = 0.18) is shown in Fig. 8 (full line) in 
comparison to the Dirac-Brueckner results from [13] (full squares) and the 
parametrizations POL6 and POL7 of the RBUU approach [9] that were found 
to optimally describe heavy-ion reactions in the energy regime up to about 1 
GeV/A. We hud the binding energy per nucleon (sa —16 MeV at pN = Po) 
to be reproduced well for a sa 0.18 which corresponds to a swelling of the 
nucleon by 18% at normal nuclear matter density. For a = 0 there is no 
minimum in F/A due to the Pauli pressure such that the swelling of the 
nucleon - which enhances the scalar attraction and reduces the vector repulsion 
- is a necessary phenomenon within the present approach to achieve proper 
binding. The resulting incompressibility K of nuclear matter amounts to about 
K sa 250 MeV. 

Since the energy per nucleon in our approach (BSD) is well in between the 
limits of POL6 and POL7, as extracted from detailed comparisons in ref. [9] for 
nucleus-nucleus collisions in the SIS energy regime, we infer that the equation 
of state generated by the model is quite realistic in the lower density (p < 3po) 
regime. Its extension to lOpo (lower part of Fig. 8), however, might still be 
questionable and has to be examined in comparison to experimental flow data 
at much higher (e.g. AGS) bombarding energies. 

The nuclear equation of state (FOS) in Fig. 8 shows no density isomer up to 
lOpo on the basis of the effective quark model adopted. The thermodynamic 


® For technical reasons we first look for the ’nucleon’ that exhibits a maximum 
quark density at a given grid point r (giving p^,Pi) and then sum up the quark 
contributions of the other ’nucleons’ (giving p^,P 2 ). The corresponding values for 
mu(r) and r°°(r) are then taken from the parametrized conhgurations displayed in 
Fig. 2. We have tested for a couple of samples that this approximate evaluation works 
quite well if ensemble averages for nuclear matter conhgurations are considered. 
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Fig. 8. Equation of state for nuclear matter; HSD (solid line), DBHF (full squares); 
RBUU results: P0L6 (dotted line), P0L7 (dashed line) from ref. [9]. 

pressure 




(30) 


furthermore, increases quadratically for pjpo > 2 and slightly levels off at high 


20 










density, but does not drop to zero in the range considered here. 

In view of the rather simple shape of the EOS in Fig. 8 and its similarity 
to the RBUU parameter sets P0L6 and P0L7 from ref. [9], it is now almost 
straight forward to ’extract’ nucleon selfenergies U§j and for the hadronic 
transport approach (1). 


2-4 Nucleon selfenergies 


The scalar and vector mean helds and in eq. (1) for nucleons now can 
be specihed along the line of ref. [7]. In order to achieve a covariant transport 
approach, which is also thermodynamically consistent [7], we parametrize the 
scalar and vector selfenergies in phase-space representation as 


UHx.p) = Ut,lx) - , 

U“{x,p) = UClx) + p ,)2 /n(^-.P')(31) 


where /n(x,p) is the nucleon phase-space distribution. In (31) the effective 
nucleon mass M*[x^p) and the kinetic momentum n^(a;,p) are given by 


M*{x,p) = + U^{x,p) , 

Il^{x,p)=p^ -U^{x,p) . (32) 

In (31) and U-lf^{x) are the local parts of the selfenergy. 


^ioc(®) = -5s 

(33) 


with 




0 4 


d‘^pfl^{x,p)f^{x,p) , 


m\ (27r)‘^ ./ 

- w{dfu,) - M*{dfu^) , 

while (Ty{{x) is obtained from the solution of 

4 


mgCTH + Ba^ J^Ca^= 


(2 


TT 


d pM*{x,p)f^{x,p). 


(34) 


(35) 
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The quasi-particle properties are defined via the mass shell constraint [7] 


h(n^n^ - M*^ 


( 36 ) 


The associated energy-momentum tensor reads 


+ {d^an{x)){d^an{x)) - {d^u^{x)){d^uf{x)) 

+ + 3-^^H + 


(jj 


H 




(27r)3 

2 

(27r)3 


(27r)3 mg 

4 gy 


Aj 

P 

A? 


(37) 




In this hadronic approach with momentum-dependent helds the Tree’ pa¬ 
rameters gg^gs^gy^gv^nis^ my , As, Ay ^ B^C allow to describe almost arbitrary 
equations of state and nucleon selfenergies. For nuclear matter at density 
the energy per nucleon is given by 


E r“ 


^ ^ - Mn , 

A pn 


(38) 


where Mn denotes the bare nucleon mass. The evaluation of T™ for 


N 


/n(p) = 20(n„)<(n' - M")0(pEi - IpI) 


(39) 


with the nucleon Fermi momentum pp then reduces to the coupled eqs. (44)- 
(49) in ref. [7] which don’t have to be repeated here. 

The key link for determining the free parameters in the hadronic model above 
now is the model independent relation for the effective quark mass as a function 
of (small) pn 


m^{pn) = mv 


1 


^ttN 

flrnl 


/5N 


(40) 


which follows from the Hellmann-Feynman theorem and the GOR relation 
(12) [50]. In (40) my is the vacuum effective quark mass from (If). In this 
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m / m 


context we show in Fig. 9 the average effective quark mass (in units of the 
vacuum mass my) (solid line) as a function of the nuclear density = p as 
obtained from the nuclear matter simulations. The effective quark mass drops 
by about 35 % at po according to (40) with XIttN = 47 MeV, and essentially 
continues with a constant slope up to about 2 X po ~ 0.33fm“^ in line with 
the Dirac-Brueckner analysis in ref. [51] (cf. also ref. [52]). The bare quark 
mass then is reached at about pN 0.6 fm It is important to note that 
the effective nucleon mass Af*(pp) (normalized to the vacuum mass) in the 
RBUU approach of ref. [9] shows the same scaling with density up to about 
Po, which is also well in line with Dirac phenomenology. 



Fig. 9. Effective mass divided by the vacuum mass as a function of the nucfeon 
density pN = p; quark mass m = m-aipN) in the HSD approach (sofid fine); nucfeon 
mass in the RBUU approach: POL6 (dotted fine), POL7 (dashed fine). 


Thus observing that the equation of state from the effective NJL model (Fig. 
8) as well as the relative scaling of the quark mass with nucleon density pN 
(Fig. 9) is very similar to the more traditional RBUU transport approach from 
refs. [7,9] at low density pN, we hx the parameters pg, ps, t/w f/v. As, Ay • • • by 
the condition 


Af*(pN,PF) _ ?7^u(Pn) 
Mn my 


(41) 
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which essentially determines the scalar selfenergy of the nucleon, as well as 
the equation of state from Fig. 8, 



7^00 

Pn 


Mn 


(42) 


up to /jn = lO/Oo- The scalar and vector nucleon selfenergies then are uniquely 
determined by eqs. (31) - (39) for arbitrary nucleon phase-space distributions. 



Fig. 10. Nucleon selfenergies Us, Uq and the Schroedinger equivalent potential Fsep 
as a function of the nucleon kinetic energy ifkin with respect to the nuclear matter 
rest frame. HSD (solid line); DBHF (full squares); exp. data from Hama et al. [53] 
(crosses). 


In Fig. 10 we compare the resulting momentum dependence of the nucleon 
selfenergies at density po with Dirac-Brueckner results from [13] (full dots). In 
the lower part of Fig. 10 the real part of the Schroedinger equivalent potential 
(SFP) 


UsEP — Us{po, P) Uo{po, P) + 

^P2 + 


2AfN 


{Us{po,Pf -Uo{po,Pf 


-\-Uo{po, P)- 


'N 


M, 


N 


(43) 
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is additionally shown (full line) in comparison to the optical potential analysis 
from Hama et ah [53] (dashed line) and Dirac-Brueckner computations from 
[13] up to momenta of 1 GeV/c. This comparison shows that the overall prop¬ 
erties of the nucleon selfenergies for Pn < 3po and Ekin < 1 GeV are reasonably 
met by our approach. 

Apart from the close analogy of our results with the a-io model at ’low’ mo¬ 
menta (cf. Fig. 10) we are especially interested in the ’high’ momentum prop¬ 
erties of the present approach, where the standard a-io model is known to fail 
signihcantly. The respective results from our present approach for the scalar 
and vector nucleon selfenergy as well as the Schroedinger equivalent optical 
potential in analogy to Fig. 10 are displayed in Fig. 11 up to relative kinetic 
energies of 15 GeV. Whereas the scalar and vector nucleon selfenergies are 
found to gradually decrease with momentum (or kinetic energy) - which is 
essentially a consequence of the cutoff Ay ~ 1.5 GeV introduced in eq. (20) - 
the Schroedinger equivalent potential exhibits a maximum of about 70 MeV at 
1 GeV and drops again for higher kinetic energy. Thus we expect the effects 
from the real part of the nucleon selfenergies to be of minor importance in 
the initial phase of nucleus-nucleus collisions at bombarding energies of a few 
GeV/A, where nucleon cascading with inelastic nucleon excitations should be 
dominant, i.e. the imaginary part of the hadron selfenergies (cf. Section 3). 



Fig. 11. Nucleon selfenergies Gs, Go and the Schroedinger equivalent potential Gsep 
as a function of the nucleon kinetic energy ifkin at normal nuclear matter density 
Po- 
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Since a transport approach for high energy nucleus-nucleus collisions also has 
to include excited states of the nucleon as well as hyperons - we include nucle¬ 
ons, A’s, N*(1440), N*(1535), A and E hyperons as well as their antiparticles - 
their respective selfenergies have to be specihed, too. As a hrst approximation 
we assume here that all baryons (made out of light (u,d) quarks) have the 
same scalar selfenergies as the nucleons; the vector selfenergy for antiparticles 
is introduced with a relative (-) sign according to time reversal® while the 
hyperons pick up a factor 2/3 according to the light quark content. 


2.5 Meson selfenergies 


Whereas the baryon selfenergies and are a necessary ingredient for a 
relativistic transport model to achieve a realistic description of hnite nuclei and 
intermediate energy nucleus-nucleus reactions, the meson selfenergies might be 
neglected in zero’th order as in conventional cascade simulations. However, in 
order to explore dynamical effects from a phase, where the chiral symmetry 
might be restored, they have to be specihed as well (on the one-loop level) e.g. 
by a suitable Lagrangian density. 

In the HSD approach, where we propagate explicitly pious, kaons, ?y’s and the 
vector mesons cu, p, and /h*(892) we assume that the pious as Goldstone 
bosons do not change their properties in the medium; we also discard selfen¬ 
ergies for the p-mesons in the present version. Thus a Lagrangian density for 
the coupled system of baryons and mesons can be written as 

£H=£B + E'C“.+4B + CB+4°B+4°t . (44) 


where Cb corresponds to the baryon Lagrangian (density) specihed in Sub¬ 
section 2.4, is the free meson Lagrangian density for a meson of type m 
and denote the meson-baryon interaction densities. The problem now is 
to hx Of^B connection with chiral symmetry restoration. 

Kaplan and Nelson [24] have shown a way how to proceed in this case. Starting 
from a SU{3)b X SU{3)r chiral Lagrangian and using chiral perturbation 
theory they write down an effective meson-baryon Lagrangian which they 
claim to be valid up to ~ 7po- Since the coefhcients in this Lagrangian are 
approximately known experimentally (within an uncertainty of about 30%), 
one can model a Lagrangian of lower complexity, but with the same properties 
on the mean-held level. Such limits lead to the following dispersion relation 
for kaons in the nuclear medium [54]: 

® This limit has to be taken with care because Teis et al. found in [3] that a sign 
change of the vector potential results in a too strong attraction for antiprotons. 
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with rri]^ denoting the bare kaon mass, /k ~ 93 MeV and Ekn ~ 350 MeV, 
while ps and p^ are the scalar and vector baryon densities, respectively. We 
thus approximate the interaction density for the K“’“-baryon system by 


^int _ 

-^K+B — 


E 


KN 


/, 


BB 


K 



+ i^(BT,.B)K+V“K+) +i^(BT‘‘B)(a,K+')K+ . (46) 

8/k 8/k 

where B, B denote the baryon spinors, and for the K“-baryon system by 


-i^(B7„B)K-*(a>‘K-) -iA(B7.-B)(a„K-*)K- (47) 

8/k 8/k 

as well as 

^ (BB) K°%° (48) 

Ik 


with the kaon holds K"*", K“ and K°. We assume the same form of /I™* for 
the K*-mesons, too. We note that Li and Ko have recently performed studies 
on the kaon dynamics at SIS energies with the same type of kaon-baryon 
interaction [55]. 

With respect to the interaction of the vector mesons p,c< 2 , with baryons we 
model according to the QCD sum rule studies by Hatsuda and Lee [56] 
as 


£" B = f — ) „y(B7»B)E-v„ 

\Po J ^ ^ 
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with Ap = sa 0.18 and A^ sa 0.025 in order to obtain a linear dependence 
of the effective meson masses on the baryon density pB = (By^B). 

With the specihcation of the meson-baryon interaction densities /1 ™b 
part of the hadron selfenergies is now fully dehned on the one-loop level. While 
the determination of the baryon selfenergies was quite involved in this Section 
in order to achieve valid approximations at low and intermediate energies as 
well as reasonable extrapolations to high baryon densities, the meson sector 
still is rather poor and will have to be improved in future. 


3 Elastic and inelastic hadron scattering 


Whereas in a fully selfconsistent relativistic transport theory the real part and 
the imaginary part of hadron selfenergies are related by means of dispersion 
relations [1,3,13], it is not justihed to employ the model selfenergies (deter¬ 
mined in Section 2) in dispersion integrals for the imaginary part because the 
inelastic scattering rate of nucleons and mesons turns out to be wrong in the 
limit of vanishing baryon density. As known from transport studies at energies 
below 2 GeV/A the elementary cross sections in eq. (1) may be approximated 
by their values in free space. Thus as a hrst step we adopt the same strategy 
and use the explicite cross-sections as in the BUU model [57] (for ^/s < 2.6 
GeV) - that have been successfully tested in the energy regime below 2 GeV/A 
bombarding energy - and by the LUND string formation and fragmentation 
model [22] (for yfs >2.6 GeV) in case of baryon-baryon collisions. For meson- 
baryon reactions we adopt a transition energy of yfs = 1.8 GeV between the 
known low energy cross sections and the LUND model. We note that the ac¬ 
tual values for the transition energies in the elementary cross sections are not 
sensitive to nucleus-nucleus collisions in the energy regime to be discussed in 
Section 4. 

In order to obtain a rough idea about the inelastic cross sections from the 
LUND string fragmentation model, we show in Fig. 12 the rapidity spectra 
for baryons, pious, kaons, p and u mesons from pp collisions at Tiab = 20 
GeV in the pp center-of-mass system. Whereas the baryons turn out to be 
located in rapidity close to the initial rapidity of the two colliding baryons, 
the meson rapidities are dominantly centered around midrapidity with a small 
contribution from the deexcitation of the baryonic constituents close to the 
incoming baryon rapidities. This general tendency has to be kept in mind 
when comparing to nucleus-nucleus collisions later on. 
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Fig. 12. Rapidity distributions for baryons, pious, kaons, r],ijj and p mesons from 
the LUND string fragmentation model for pp collisions at Tiab = 20 GeV. 

The implementation of the LUND string formation and fragmentation model 
[22] - which describes the free transition probabilities - in a covariant transport 
theory implies to use a time scale to transform the cross-sections to collision 
rates and particle production rates. An appropriate time scale is given by a 
string formation time Up, which denotes the time between the formation and 
fragmentation of the string in the individual hadron-hadron center-of-mass 
system for a particle of rapidity j/cm = 0. Due to covariance this time should 
be also related to the spatial extension of the interacting hadrons which on 
average gives Up ~ 0.8 fm/c. The sensitivity of the proton rapidity spectra 
dNIdy to the actual value of Up ist shown in Fig. 13 for central collisions of 
“^^Ca + “^^Ca at 30 GeV/A. It is seen that Up controls essentially the rapidity 
distribution at midrapidity and at projectile and target rapidity, i.e. the baryon 
stopping in relativistic nucleus-nucleus collisions. We will adopt rp=0.8 fm/c 
for the calculations to be presented in Section 4; similar values are also used 
in the RQMD approach [58]. 
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Fig. 13. The proton rapidity distribution for central collisions of “^^Ca + “^^Ca at 30 
GeV/A using rF= 0.5, 0.6, 0.8, 1.0 fm/c, respectively. 

In view of the ’chiral’ dynamics addressed in this work, however, especially 
the productions rates of mesons should change at high baryon density [59] 
due to the reduced masses involved. Unfortunately, the actual meson scalar 
and vector selfenergies are quite a matter of debate and depend on the model 
parameters of the Lagrangians employed. Following the approach by Kaplan 
and Nelson [24] (c.f. Section 2.5), the average mass of a K"*" K“ pair is expected 
to follow approximately 

mn+K- ~ rny{l — O.lQp^/po) (50) 


because the vector interaction drops out due to contributions with opposite 
sign and the term in is rather small. Furthermore, according to the meson- 
baryon interaction densities (49) the in-medium mass for p, io and $ mesons 
- following Hatsuda and Lee [56]- can be approximated by 

m(pN) ~ ?7^v(l - Apn//Oo) (51) 


with A sa 0.18 for p, lo and A sa 0.025 for $ mesons. The weak dependence of 
the $ meson mass here is a consequence of the weak coupling of the strange 
quark to the light (u,d) quarks which dominantly make up the baryon density. 

In view of the substantial uncertainties of the meson selfenergies especially 
at high density we here use this more pragmatic model which does not claim 
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fundamental evidence ^. Whereas the pion as a Goldstone boson is assumed 
not to change substantially with baryon density and temperature in the energy 
regime addressed, the kaons, K*’s, p’s, cu’s and ^’s are assumed to change their 
masses as displayed in Fig. 14 roughly in line with ref. [24,56] as pointed out 
above. The hnal values achieved at high baryon density are determined by the 
bare quark mass content of the mesons. 



Fig. 14. Parametrization of the effective meson masses - normalized to the vacuum 
masses niy - versus the baryon density used in the extended string fragmentation 
model (HSD). 


As an example for the effects to be expected at high baryon density we show in 
Fig. 15 the rapidity spectra of kaons and p’s for pp collisions at Tiab =20 GeV 
from the string fragmentation model that incorporates the density dependent 
meson masses from Fig. 14® . It is clearly seen that a dropping of the meson 
masses leads to a substantial enhancement of the K“’“K“ and p yields and to a 
widening of their rapidity distribution in the individual center-of-mass system. 
If such phenomena can be seen in comparison to experimental data will be 
investigated in the next Section. 


^ These assumptions about meson properties at high baryon density can only be 
controlled in confrontation with sensitive experimental data. 

® The total four-momentum is conserved in the ’chiral’ string fragmentation model. 
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Fig. 15. Differential K+ + K and p rapidity distributions from the ’chiral’ string 
fragmentation model at 1.0, 0.7 and 0.4 X the bare masses for pp collisions at 
Fiab=20 GeV. 

4 Heavy-ion collisions 


The relativistic transport approach (HSD) outlined in Sections 2 and 3 now 
is applied to nucleus-nucleus collisions from the SIS to the SPS energy regime 
with particular emphasis on rapidity distributions and particle spectra to con¬ 
trol the stopping achieved in these reactions. The explicit numerical implemen¬ 
tation of the selfenergies and collisions rates is performed in close analogy to 
[9,57,60,61] and does not have to be repeated here. We note that the total 
conservation of energy and momentum throughout the time evolution is con¬ 
served on the 2 % level for central An + An collisions and even better for 
peripheral or light-ion induced reactions. 

As a hrst example we show in Fig. 16 the transverse 7r°-spectra from Ar + 
Ca collisions at 1.5 GeV/A in comparison to the data of Berg et ah [62] as 
a characteristic system at SIS energies. Since at these energies the present 
approach is close to the results achieved with the former BUU model [63], the 
reproduction of the data is of similar quality. We thus conclude that the ’low 
energy dynamics’ involving essentially nucleons, A’s and pious is reasonably 
well included in our transport calculations. 


^.1 Stopping in high-energy nucleus-nucleus collisions 


The next system addressed is Si + A1 at 14.6 GeV/A, i.e. the AGS energy 
regime. The computed rapidity distribution of protons and tt“’“- mesons for 
b = 1.5 fm is compared in Fig. 17 to the data from ref. [64]. Whereas the 
proton rapidity distribution turns out to be quite flat in rapidity y due to 
proton rescattering, the pion rapidity distribution is essentially of gaussian 
shape which reflects the pion rapidity spectrum from the string fragmentation 
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Fig. 16. Calculated transverse 7r°-spectra for Ar + Ca at 1.5 GeV/A (full line) in 
comparison to the experimental data from ref. [62]. 


model outlined in Section 3 (cf. Fig. 12). Similar to SIS energies [1,2] the 
proton rapidity distribution is insensitive to variations of the nucleon scalar 
and vector mean helds within the numerical accuracy. 

In analogy to Fig. 16 we show in Fig. 18 the calculated transverse mass-spectra 
of TT“’“-mesons for Si + A1 at 14.6 GeV/A (solid lines) in comparison to the 
experimental data from ref. [64]. The overall agreement for lab. rapidities of 
y = 0.9, 1.7 and 2.7 seems to indicate that the general reaction dynamics is 
well reproduced within the HSD approach. 

The flat proton rapidity spectrum in Fig. 17 might lead to the interpretation 
that there is a substantial amount of stopping in the light system Si + Ah 
This, however, has to be taken with care because the actual snapshots of the 
baryon density distribution from our computations shown in Fig. 19 (l.h.s.) 
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ycm 

Fig. 17. Calculated proton and 7r+ rapidity distribution (liistrograms) for a central 
14.6 GeV/A Si+Al collision in comparison to the data from [64] (full dots). 

as well the phase-space distribution (r.h.s.) 

f{z,p,;t) = {2Tr)-^J2 / dr^dp^fb{r^,z,p^,p,-t) , (52) 

where denotes a sum over all baryon species, indicate a dominant trans¬ 
parency for the light system. This is essentially due to the large surface of the 
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Fig. 18. Comparison of the calculated transverse mass spectra of vr■’“-mesons for Si 
-|- A1 at 14.6 GeV/A with the experimental data from ref. [64] for rapidities y=0.9, 
1.7, 2.7 in the laboratory system. 

two light nuclei with a nucleon-nucleon collision probability less than 1. Fur¬ 
thermore, the time evolution in momentum space (middle column) shows that 
the system is far from kinetic equilibrium in the baryon degrees of freedom in 
the hnal state. 

The amount of stopping at AGS energies is more clearly pronounced for cen¬ 
tral An + An reactions as displayed in Fig. 20 for the proton and tt~ rapidity 
distributions. Though the pion rapidity spectrum does not differ very much 
in shape from that of the Si + A1 system in Fig. 17 at first sight, the time 
evolution of the baryon distribution in coordinate space, momentum space 
and phase space (Fig. 21) for Au+Au at 14.6 GeV/A shows a clear approach 
versus equilibration. However, the coordinate space evolution indicates a dom¬ 
inant longitudinal expansion which is also reflected in the baryon momentum 
distribution that does not show full isotropy. Detailed experimental data and 
related comparisons, however, will become available soon at the energy of 10.8 
GeV/A [65]. We note that the proton rapidity spectrum for central An + An 
collisions at this energy shows a similar amount of stopping as the RQMD 
approach [58]. 

We continue our comparison to experimental data with the system S + S at 200 
GeV/A, i.e. the SPS regime (Fig. 22). Though the experimental proton and tt~ 
rapidity spectra (from [66]) are approximately reproduced, we cannot conclude 
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Fig. 19. Baryon density distribution (left column), momentum space (middle col¬ 
umn) and phase-space distribution (right column) for a 14.6 GeV/A Si+Al collision 
at b = 1 fm for various times in fm/c. 
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Fig. 20. Proton (dashed line) and 7r“ rapidity distribution (full line) for a central 
10.8 GeV/A An + An collision. 

on the general applicability of our approach at SPS energies because also more 
simple models like HIJING or VENUS - with a less amount of rescattering - 
can reproduce the data in a similar way [67]. A way out of this problem is to 
analyze the system Pb + Pb at 153 GeV/A (Fig. 23) that has recently been 
studied experimentally at the SPS. Our computed proton rapidity spectrum 
for central collisions shows no dip at midrapidity as in HIJING or VENUS 
simulations [67] but a flat spectrum similar to RQMD simulations [58]. On 
the other hand, the pion rapidity distributions are very similar to the S + S 
case, however, enhanced by about a factor of 7.5. 


4.2 Probing chiral symmetry restoration 


The problem of chiral symmetry restoration can be investigated e.g. via the 
K“’“/7r“’“ ratio as addressed in [59] since the kaon production should be en¬ 
hanced due to the reduced in-medium mass (cf. Fig. 15). For this purpose we 
show in Table 1 the calculated K“’“/7r“’“ yields for the systems p + p. Si + Al, 
Si + An at 14.6 GeV/A and An + An in comparison with the experimental 
data for two different szenarios. The first column represents the results of a 
simulation where only the bare masses of the mesons have been considered 
in the string fragmentation approach (HSD) and all mesons are propagated 
as free particles whereas the second column results from density-dependent 
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Fig. 22. Proton and vr rapidity distribution for a central 200 GeV/A S + S collision 
in comparison to the data of ref. [66]. 









exp.ratio 

without kaon 

selfenergies 

with kaon 

selfenergies 

p + p 

0.08 

0.08 

0.08 

Si + Al 

0.13 

0.09 

0.12 

Si + Cu 

0.16 

0.1 

0.15 

Si -|- Au 

0.19 

0.11 

0.16 

Au -K Au 

0.22 

0.12 

0.21 


Table 1. The K+/7r+ ratio for p + p, Si + Al, Si + Au and An + An collisions at 
14.6 GeV/A in comparison to the data from ref. [68]. 


meson masses as described in Section 3. According to their drop in mass 
Amh(r,f) = 0 mesons are propagated in their time-dependent 

scalar mean held which couples linearly to the baryon density. Thus 

their momenta are decreased dynamically during the expansion of the hadronic 
system and the energy to become asymptotically on-shell is extracted from 
the collective motion. It is clearly seen that for density-dependent K"*" masses 
the ratio is strongly enhanced for the heavier systems as seen experimentally. 
However, this enhancement could also be attributed to a closer approach to 
chemical equilibrium as advocated in ref. [68] which might be achieved due to 
enhanced hadronic cross sections in the dense medium. Whereas in principle 
the coupled transport equations (1) also describe the approach towards chem¬ 
ical equilibrium for large systems, it is not yet clear if all the proper reactions 
rates are presently included in our simulations such that no hnal evidence on 
chiral symmetry restoration can be extracted so far from the K“’“/7r“’“ ratio. 


The medium modihcations of the p-meson are most efficiently probed by dilep¬ 
ton spectroscopy [57,59,63] since due to its short lifetime the p-meson has a 
good chance to decay in the dense baryonic environment. According to Fig. 15 
we expect a substantial enhancement of dileptons in the invariant mass range 
0.4 GeV < M < 0.7 GeV in nucleus-nucleus collisions as compared to p + 
A collisions due to a shift in the p-mass spectrum and an enhanced p-meson 
production in the dense medium especially via 7r“’“7r“ annihilation [69,70]. In 
fact, hrst computations for dilepton production within the HSD approach [70] 
show that the enhancement of dileptons in central S + Au collisions at 200 
GeV/A (reported by the CERES-collaboration [71]) might be explained by 
the ’chiral’ dynamics proposed in Sections 2 and 3. 
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5 Summary 


In this work we have presented a relativistic transport approach for hadrons 
(denoted by HSD® ) where the underlying (real parts of) nucleon selfenergies 
have been determined on the basis of an effective NJL-type Lagrangian for the 
quark degrees of freedom with a chiral invariant interaction density. Starting 
with a local color current interaction we have developed a model for spin 
and isospin averaged color neutral states on the basis of the experimental 
electromagnetic formfactor of the proton. The parameters in our model, which 
are all hxed by physical quantities are Gs and As for the scalar part, Gy and 
Ay for the vector part and a, which describes the swelling of the nucleon in 
nuclear matter. The physical quantities which are sufficiently well met are: 
the averaged nucleon mass Mn, the scalar vacuum condensate the 

pion-nucleon S-term, the nuclear equation of state (minimum at p = po with 
-16 MeV binding energy) and the Schroedinger equivalent potential ffisEP for 
nucleons. Due to the scalar-vector nature of the quark-quark interaction the 
nucleon selfenergies (htted to the NJL results) are also close to those of the 
cr-io model of Walecka [40] in the low momentum and low density regime. 

Whereas the real part of the nucleon selfenergies has been determined from a 
quark oriented model to allow for reasonable extrapolations to the high density 
regime, the kaon, p, io and ^ meson selfenergies are hxed in line with more 
simple Lagrangian densities (Section 2.5). In this respect the meson sector 
will need further improvement in future. However, especially at bombarding 
energies as high as 200 GeV/A, the imaginary part of the hadron selfener¬ 
gies is more important. In the HSD approach the respective transition rates 
have been adopted from the LUND string fragmentation model [22] where the 
meson masses (except the pion and rj) have been scaled in line with chirally 
invariant interaction densities. This more pragmatic model, of cause, has less 
founded reliability and thus one has to justify its applicability or inadequacy 
in comparison to experimental data. 

As a hrst step in this direction we have applied our relativistic transport ap¬ 
proach to nucleus-nucleus collisions from SIS to SPS energies. Whereas the 
proton and pion rapidity distributions and transverse pion spectra look rea¬ 
sonably well for the systems studied experimentally, a clear signature for the 
chiral symmetry restoration could not unambiguously be established so far. 
This is because the strangeness enhancement observed experimentally at AGS 
and SPS energies might also be due to chemical equilibration or e.g. color-rope 
formation [58]. A better probe should be provided by dilepton spectroscopy 
in the invariant mass regime from 0.4 - 0.8 GeV [57] since the p-meson pre¬ 
dominantly decays in the dense medium. In fact, hrst computations on e“’“e“ 
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production at SPS energies suggest that the dilepton enhancement seen by 
the CERES collaboration might be due to a dropping p-mass in the medium 

[70]. 

The nuclear equation of state computed within our approach (Eig. 8) shows 
no density isomer up to p sa lOpo- This prediction is essentially due to the 
fact that scalar and vector baryon selfenergies are approximately of the same 
order of magnitude, but different sign, up to p sa 4po and the repulsive vector 
interaction takes over at even higher p together with the kinetic energy per 
nucleon. A density isomer thus can only occur if the vector coupling itself 
decreases at high baryon density or temperature. Some arguments in this di¬ 
rection have recently been proposed by Brown and Rho [72] and investigated 
in a model study by Li and Ko [73]. A clarihcation of this problem e.g. should 
be achieved by experimental data on the baryon flow as a function of pro¬ 
jectile/target mass and bombarding energy in the energy regime in between 
SIS and SPS thus allowing for a closer look at the EOS at Very high’ baryon 
density. 


The authors acknowledge valuable and inspiring discussions througout this 
work with C. M. Ko, U. Mosel, H. Stocker, S. Teis and Gy. Wolf. They are 
also grateful to T. Maruyama for an earlier version of the relativistic transport 
code used in the analysis of nucleus-nucleus collisions in the energy regime 
below 1 GeV/A. 
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